library(foreign);library(lme4);library(arm);library(xtable);library(car);library(ggplot2);library(stargazer);library(survey)

model <- lm(scsupx~IDEODIF2+SCJOB2+rlawx+tolx+orderx+SCKNOWX5+partyid+age+hispanic+black+degree+gender+ownhomed+attend+BAGAIN2,data);summary(model)

res <- read.csv("results.csv")
###HELD AT MEANS
#BASELINE
fv <- res$coef[1]*res$mean[1] + res$coef[2]*res$mean[2] + res$coef[3]*res$mean[3] + res$coef[4]*res$mean[4] + res$coef[5]*res$mean[5] + res$coef[6]*res$mean[6] + res$coef[7]*res$mean[7] + res$coef[8]*res$mean[8] + res$coef[9]*res$mean[9] + res$coef[10]*res$mean[10] + res$coef[11]*res$mean[11] + res$coef[12]*res$mean[12] + res$coef[13]*res$mean[13] + res$coef[14]*res$mean[14] + res$coef[15]*res$mean[15] + res$coef[16]

##IQR
fvMEANlow1 <- res$coef[1]*.0177777777777778 + res$coef[2]*res$mean[2] + res$coef[3]*res$mean[3] + res$coef[4]*res$mean[4] + res$coef[5]*res$mean[5] + res$coef[6]*res$mean[6] + res$coef[7]*res$mean[7] + res$coef[8]*res$mean[8] + res$coef[9]*res$mean[9] + res$coef[10]*res$mean[10] + res$coef[11]*res$mean[11] + res$coef[12]*res$mean[12] + res$coef[13]*res$mean[13] + res$coef[14]*res$mean[14] + res$coef[15]*res$mean[15] + res$coef[16]

fvMEANhigh1 <- res$coef[1]*.4444444444444444 + res$coef[2]*res$mean[2] + res$coef[3]*res$mean[3] + res$coef[4]*res$mean[4] + res$coef[5]*res$mean[5] + res$coef[6]*res$mean[6] + res$coef[7]*res$mean[7] + res$coef[8]*res$mean[8] + res$coef[9]*res$mean[9] + res$coef[10]*res$mean[10] + res$coef[11]*res$mean[11] + res$coef[12]*res$mean[12] + res$coef[13]*res$mean[13] + res$coef[14]*res$mean[14] + res$coef[15]*res$mean[15] + res$coef[16]

fvMEANhigh1-fvMEANlow1

fvMEANlow2 <- res$coef[1]*res$mean[1] + res$coef[2]*.3333333333333333 + res$coef[3]*res$mean[3] + res$coef[4]*res$mean[4] + res$coef[5]*res$mean[5] + res$coef[6]*res$mean[6] + res$coef[7]*res$mean[7] + res$coef[8]*res$mean[8] + res$coef[9]*res$mean[9] + res$coef[10]*res$mean[10] + res$coef[11]*res$mean[11] + res$coef[12]*res$mean[12] + res$coef[13]*res$mean[13] + res$coef[14]*res$mean[14] + res$coef[15]*res$mean[15] + res$coef[16]

fvMEANhigh2 <- res$coef[1]*res$mean[1] + res$coef[2]*.6666666666666666 + res$coef[3]*res$mean[3] + res$coef[4]*res$mean[4] + res$coef[5]*res$mean[5] + res$coef[6]*res$mean[6] + res$coef[7]*res$mean[7] + res$coef[8]*res$mean[8] + res$coef[9]*res$mean[9] + res$coef[10]*res$mean[10] + res$coef[11]*res$mean[11] + res$coef[12]*res$mean[12] + res$coef[13]*res$mean[13] + res$coef[14]*res$mean[14] + res$coef[15]*res$mean[15] + res$coef[16]

fvMEANhigh2-fvMEANlow2


##1.96 SDS 
fvMEANlow1 <- res$coef[1]*(.2238022-1.96*.2411775) + res$coef[2]*res$mean[2] + res$coef[3]*res$mean[3] + res$coef[4]*res$mean[4] + res$coef[5]*res$mean[5] + res$coef[6]*res$mean[6] + res$coef[7]*res$mean[7] + res$coef[8]*res$mean[8] + res$coef[9]*res$mean[9] + res$coef[10]*res$mean[10] + res$coef[11]*res$mean[11] + res$coef[12]*res$mean[12] + res$coef[13]*res$mean[13] + res$coef[14]*res$mean[14] + res$coef[15]*res$mean[15] + res$coef[16]

fvMEANhigh1 <- res$coef[1]*(.2238022+1.96*.2411775) + res$coef[2]*res$mean[2] + res$coef[3]*res$mean[3] + res$coef[4]*res$mean[4] + res$coef[5]*res$mean[5] + res$coef[6]*res$mean[6] + res$coef[7]*res$mean[7] + res$coef[8]*res$mean[8] + res$coef[9]*res$mean[9] + res$coef[10]*res$mean[10] + res$coef[11]*res$mean[11] + res$coef[12]*res$mean[12] + res$coef[13]*res$mean[13] + res$coef[14]*res$mean[14] + res$coef[15]*res$mean[15] + res$coef[16]

fvMEANlow1-fvMEANhigh1

fvMEANlow2 <- res$coef[1]*res$mean[1] + res$coef[2]*(.5757985-1.96*.2220815) + res$coef[3]*res$mean[3] + res$coef[4]*res$mean[4] + res$coef[5]*res$mean[5] + res$coef[6]*res$mean[6] + res$coef[7]*res$mean[7] + res$coef[8]*res$mean[8] + res$coef[9]*res$mean[9] + res$coef[10]*res$mean[10] + res$coef[11]*res$mean[11] + res$coef[12]*res$mean[12] + res$coef[13]*res$mean[13] + res$coef[14]*res$mean[14] + res$coef[15]*res$mean[15] + res$coef[16]

fvMEANhigh2 <- res$coef[1]*res$mean[1] + res$coef[2]*(.5757985+1.96*.2220815) + res$coef[3]*res$mean[3] + res$coef[4]*res$mean[4] + res$coef[5]*res$mean[5] + res$coef[6]*res$mean[6] + res$coef[7]*res$mean[7] + res$coef[8]*res$mean[8] + res$coef[9]*res$mean[9] + res$coef[10]*res$mean[10] + res$coef[11]*res$mean[11] + res$coef[12]*res$mean[12] + res$coef[13]*res$mean[13] + res$coef[14]*res$mean[14] + res$coef[15]*res$mean[15] + res$coef[16]

fvMEANlow2-fvMEANhigh2

###At each IQR
#0
res$coef[1]*0 + res$coef[2]*res$mean[2] + res$coef[3]*res$mean[3] + res$coef[4]*res$mean[4] + res$coef[5]*res$mean[5] + res$coef[6]*res$mean[6] + res$coef[7]*res$mean[7] + res$coef[8]*res$mean[8] + res$coef[9]*res$mean[9] + res$coef[10]*res$mean[10] + res$coef[11]*res$mean[11] + res$coef[12]*res$mean[12] + res$coef[13]*res$mean[13] + res$coef[14]*res$mean[14] + res$coef[15]*res$mean[15] + res$coef[16]

#.25 & .33
res$coef[1]*.0177777777777778 + res$coef[2]*res$mean[2] + res$coef[3]*res$mean[3] + res$coef[4]*res$mean[4] + res$coef[5]*res$mean[5] + res$coef[6]*res$mean[6] + res$coef[7]*res$mean[7] + res$coef[8]*res$mean[8] + res$coef[9]*res$mean[9] + res$coef[10]*res$mean[10] + res$coef[11]*res$mean[11] + res$coef[12]*res$mean[12] + res$coef[13]*res$mean[13] + res$coef[14]*res$mean[14] + res$coef[15]*res$mean[15] + res$coef[16]


#.5
res$coef[1]*.16 + res$coef[2]*res$mean[2] + res$coef[3]*res$mean[3] + res$coef[4]*res$mean[4] + res$coef[5]*res$mean[5] + res$coef[6]*res$mean[6] + res$coef[7]*res$mean[7] + res$coef[8]*res$mean[8] + res$coef[9]*res$mean[9] + res$coef[10]*res$mean[10] + res$coef[11]*res$mean[11] + res$coef[12]*res$mean[12] + res$coef[13]*res$mean[13] + res$coef[14]*res$mean[14] + res$coef[15]*res$mean[15] + res$coef[16]

#.67
res$coef[1]*.284444444444444 + res$coef[2]*res$mean[2] + res$coef[3]*res$mean[3] + res$coef[4]*res$mean[4] + res$coef[5]*res$mean[5] + res$coef[6]*res$mean[6] + res$coef[7]*res$mean[7] + res$coef[8]*res$mean[8] + res$coef[9]*res$mean[9] + res$coef[10]*res$mean[10] + res$coef[11]*res$mean[11] + res$coef[12]*res$mean[12] + res$coef[13]*res$mean[13] + res$coef[14]*res$mean[14] + res$coef[15]*res$mean[15] + res$coef[16]

#.75
res$coef[1]*.4444444 + res$coef[2]*res$mean[2] + res$coef[3]*res$mean[3] + res$coef[4]*res$mean[4] + res$coef[5]*res$mean[5] + res$coef[6]*res$mean[6] + res$coef[7]*res$mean[7] + res$coef[8]*res$mean[8] + res$coef[9]*res$mean[9] + res$coef[10]*res$mean[10] + res$coef[11]*res$mean[11] + res$coef[12]*res$mean[12] + res$coef[13]*res$mean[13] + res$coef[14]*res$mean[14] + res$coef[15]*res$mean[15] + res$coef[16]

#1
res$coef[1]*1 + res$coef[2]*res$mean[2] + res$coef[3]*res$mean[3] + res$coef[4]*res$mean[4] + res$coef[5]*res$mean[5] + res$coef[6]*res$mean[6] + res$coef[7]*res$mean[7] + res$coef[8]*res$mean[8] + res$coef[9]*res$mean[9] + res$coef[10]*res$mean[10] + res$coef[11]*res$mean[11] + res$coef[12]*res$mean[12] + res$coef[13]*res$mean[13] + res$coef[14]*res$mean[14] + res$coef[15]*res$mean[15] + res$coef[16]



###HELD AT MEDIANS
#BASELINE
fv <- res$coef[1]*res$median[1] + res$coef[2]*res$median[2] + res$coef[3]*res$median[3] + res$coef[4]*res$median[4] + res$coef[5]*res$median[5] + res$coef[6]*res$median[6] + res$coef[7]*res$median[7] + res$coef[8]*res$median[8] + res$coef[9]*res$median[9] + res$coef[10]*res$median[10] + res$coef[11]*res$median[11] + res$coef[12]*res$median[12] + res$coef[13]*res$median[13] + res$coef[14]*res$median[14] + res$coef[15]*res$median[15] + res$coef[16]

fvMEDIANlow1 <- res$coef[1]*.0177777777777778 + res$coef[2]*res$median[2] + res$coef[3]*res$median[3] + res$coef[4]*res$median[4] + res$coef[5]*res$median[5] + res$coef[6]*res$median[6] + res$coef[7]*res$median[7] + res$coef[8]*res$median[8] + res$coef[9]*res$median[9] + res$coef[10]*res$median[10] + res$coef[11]*res$median[11] + res$coef[12]*res$median[12] + res$coef[13]*res$median[13] + res$coef[14]*res$median[14] + res$coef[15]*res$median[15] + res$coef[16]
fvMEDIANhigh1 <- res$coef[1]*.4444444444444444 + res$coef[2]*res$median[2] + res$coef[3]*res$median[3] + res$coef[4]*res$median[4] + res$coef[5]*res$median[5] + res$coef[6]*res$median[6] + res$coef[7]*res$median[7] + res$coef[8]*res$median[8] + res$coef[9]*res$median[9] + res$coef[10]*res$median[10] + res$coef[11]*res$median[11] + res$coef[12]*res$median[12] + res$coef[13]*res$median[13] + res$coef[14]*res$median[14] + res$coef[15]*res$median[15] + res$coef[16]

fvMEDIANhigh1-fvMEDIANlow1

fvMEDIANlow2 <- res$coef[1]*res$median[1] + res$coef[2]*.3333333333333333 + res$coef[3]*res$median[3] + res$coef[4]*res$median[4] + res$coef[5]*res$median[5] + res$coef[6]*res$median[6] + res$coef[7]*res$median[7] + res$coef[8]*res$median[8] + res$coef[9]*res$median[9] + res$coef[10]*res$median[10] + res$coef[11]*res$median[11] + res$coef[12]*res$median[12] + res$coef[13]*res$median[13] + res$coef[14]*res$median[14] + res$coef[15]*res$median[15] + res$coef[16]
fvMEDIANhigh2 <- res$coef[1]*res$median[1] + res$coef[2]*.6666666666666666 + res$coef[3]*res$median[3] + res$coef[4]*res$median[4] + res$coef[5]*res$median[5] + res$coef[6]*res$median[6] + res$coef[7]*res$median[7] + res$coef[8]*res$median[8] + res$coef[9]*res$median[9] + res$coef[10]*res$median[10] + res$coef[11]*res$median[11] + res$coef[12]*res$median[12] + res$coef[13]*res$median[13] + res$coef[14]*res$median[14] + res$coef[15]*res$median[15] + res$coef[16]

fvMEDIANhigh2-fvMEDIANlow2





####Plot Sophistication Interaction
ruler <- seq(0,1,length.out=100)
meff <- -.0319051 + ruler*-.0810474
se <- sqrt(.0590412^2 + ruler*ruler*.0872809*.0872809+2*ruler*-.00462096)

par(mfrow=c(1,1))
plot(ruler,meff,type="l",ylim=c(-.25,.1),xlab="Liberal                                                                Knowledge                                                            Conservative",ylab="Marginal Effect",main="Marginal Effect of Ideological Disagreement on Legitimacy, \n by Supreme Court Knowledge",lwd=2)
lines(ruler,meff+1.96*se,lty=2)
lines(ruler,meff-1.96*se,lty=2)
abline(h=0,col="red",lwd=2)

####Plot Effect of ideology across supreme court
ruler <- seq(0,1,length.out=100)
meff <- -.0932389  + ruler*.112127 
se <- sqrt(.0344686^2 + ruler*ruler*.0725708*.0725708+2*ruler*-.00200909)

plot(ruler,meff,type="l",ylim=c(-.2,.2),xlab="Liberal                                             Respondent Ideological Placement                                        Conservative",ylab="Marginal Effect",main="Marginal Effect of Supreme Court Ideological Placement on Legitimacy, \n by Respondent Ideological Placement",lwd=2)
lines(ruler,meff+1.96*se,lty=2)
lines(ruler,meff-1.96*se,lty=2)
abline(h=0,col="red",lwd=2)

####Plot Effect of supreme court across ideology
ruler <- seq(0,1,length.out=100)
meff <- -.0291827 + ruler*.112127 
se <- sqrt(.0398192^2 + ruler*ruler*.0725708*.0725708+2*ruler*-.00216774 )

plot(ruler,meff,type="l",ylim=c(-.5,.5),xlab="Supreme Court Ideological Placement",ylab="Marginal Effect",main="Marginal Effect of Respondent Ideological Placement on Legitimacy, \n by Supreme Court Ideological Placement",lwd=2)
lines(ruler,meff+1.96*se,lty=2)
lines(ruler,meff-1.96*se,lty=2)
abline(h=0,col="red",lwd=2)



